Statistics of thermal to shot noise crossover in 
chaotic cavities 



Vladimir Al. Osipov^ and Eugene Kanzieper^ 

^ Fachbereich Physik, Universitat Duisburg-Essen, Duisburg D-47057, Germany 
Department of Applied Mathematics, H.I.T. - Holon Institute of Technology, 
Holon 58102, Israel 

E-mail: vladimir . oslpov{Buni-due . de and eugene.kaiizieperahit.ac.il 

Abstract. Recently formulated integrable theory of quantum transport [Osipov 
and Kanzieper, Phys. Rev. Lett. 101, 176804 (2008)] is extended to describe 
sample-to-sample fluctuations of the noise power in chaotic cavities with broken 
time-reversal symmetry. Concentrating on the universal transport regime, we de- 
termine dependence of the noise power cumulants on the temperature, applied bias 
voltage, and the number of propagating modes in the leads. Intrinsic connection 
between statistics of thermal to shot noise crossover and statistics of Landauer 
conductance is revealed and briefly discussed. 
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1. Introduction: Thermal versus shot noise 

The charge transfer through a phase-coherent cavity exhibiting chaotic classical 
dynamics is a random process influenced by discreteness of the electron charge e 
and the quantum nature of electrons (Blantcr and Biittikcr 2000, Imry 2002, Martin 
2005). Fluctuations of charge transmitted during a fixed time interval or, equivalently, 
fluctuations SI{t) of current around its mean are quantified by the noise power 

/ + 00 
dt {SI{t + to)SI{to))t^, (1.1) 
-oo 

where the brackets (• • ■)^^^ indicate the time averaging. 

At temperatures 6 = ksT which are much larger than a bias voltage v = 
eV applied to the cavity {9 » v), the current fluctuations are dominated by 
the equilibrium thermal noise, also known as Johnson-Nyquist noise. Caused by 
fluctuating occupation numbers in a flow of carriers injected into cavity from electronic 
reservoirs, thermal noise extends over all frequencies up to the quantum limit 6/h. In 
the absence of electron-electron interactions, its power at zero bias voltage {v = 0) is 
related to the scattering matrix § of the system composed of the cavity and the leads 
(Khlus 1987, Lesovik 1989, Biittiker 1990, Martin and Landauer 1992, Biittiker 1992): 

TthW =4eGotr(eiSe2St). (1.2) 

Here, Go = e'^/h is the conductance quantum. The projection matrices Ci,2 encode 
the information about particular cavity-lcad geometry and will be specified later on. 

In the opposite limit of low temperatures {6 <C v), the current fluctuations are 
still significant even though the flow of incident electrons is essentially noiseless. In 
this temperature regime, noncquilibrium current fluctuations (known as a shot noise) 
exist because of (i) the granularity of the electron charge e and (ii) the stochastic 
nature of electron scattering inside the cavity which splits the electron wave into two 
or more; partial waves leaving the cavity through diffcrcnit exits. It is this "uncertainty 
of not knowing where the electron came from and where it will go to" (Oberholzer 
et al 2002) that makes the transmitted charge to fluctuate. At zero temperature, the 
scattering matrix approach brings the shot noise power in the form 

3'shot(i^) = 2t; Go [tr(eiSe2§t) - tr(eiSe2St)2] . (1.3) 

At finite temperatures, both sources of noise are operative, the total noise CP(^, v) 
being a complicated function of temperature and bias voltage p. 

9(9, v) = A9 Go (tr(eiSe2§t)2 + ^ coth (^) [tr(ei§e2St) - tr(eiSe2St)2] ) . (1.4) 

Equation (1.4) suggests that the crossover from thermal noise 5'th(^^) = "^{9, 0) to shot 
noise J'shot(^^) = ^(0,^) depends in a sensitive way on scattering properties of the 
cavity and the leads incorporated in the scattering matrix S. Since chaotic scattering 
of electrons inside the cavity induces fluctuations of S-matrix (Bliimel and Smilansky 
1990), the noise power 'S'{9,v) fiuctuates, too. 

So far, the thermal to shot noise crossover has only been studied at the level of 
average noise power. For the two-terminal scattering geometry comprised of the cavity 
attached to outside reservoirs (kept at temperature 9) via two leads supporting TVl 

X Equation (1.4) disregcirds the low-frequency 1// noise that can efficiently be filtered out in 
experiments. 
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and Nji propagating modes, respectively, the average noise power equals (Blanter and 
Sukhorukov 2000, Oberholzer et al 2001, Savin and Sommers 2006) 



(T(^,i;))g = (Tth)s 

where 



(1.5) 



is the average equilibrium thermal noise power, and the thermodynamic function 

/^=/3coth/3-l (1.7) 

depends on the ratio /3 = v/29 between the bias voltage v and the temperature 9. 
Equations (1.5) and (1.6) hold for cavities with broken time reversal symmetry §. 
Derived for the universal transport regime (Beenakker 1997, Richter and Sieber 2002, 
Miiller et al 2007) emerging in the hmit td » m (Agam et al 2000), where td is 
the average electron dwell time and te is the Ehrenfest time (the time scale where 
quantum effects set in), the above prediction has been confirmed in a remarkable series 
of experiments (Oberholzer et al 2001, Cronei al 2001, Oberholzer et al 2002). 

In this paper, we examine statistics of the thermal to shot noise crossover. The 
latter, contained in the distribution function of the noise power CP(0, v) or. cquivalcntly, 
in its cumulants ((J'^)), can effectively be described within the framework of integrable 
theory of quantum transport formulated by Osipov and Kanzieper (2008). Let us stress 
that recent experimental studies (Flindt et al 2009) of quantum noise fluctuations 
in nanoscale conductors (which concentrated on detection of higher cumulants of 
noise) suggest that testing our predictions may be feasible within the current limits 
of nanotechnology. 



2. Integrable theory of noise power fluctuations 

In what follows, we consider chaotic cavities with broken time-reversal symmetry 
which are probed, via ballistic point contacts, by two (left and right) leads; the 
leads supporting A^l and propagating modes, respectively, are further coupled to 
external reservoirs kept at the temperature 9. This scattering geometry corresponds 
to the projection matrices Ci 2 of the form 

see Eqs. (1.2) - (1.4). 

2.1. Joint cumulants of Landauer conductance and noise power 

The starting point of our analysis is the joint cumulant generating function (JCGF) 

Jn{z, w) = {exp{-z G/Go) exp{-w secVEiN) (2-2) 

of the Landauer conductance G = Gotr(eiSC2§^) and the noise power 7(9, v) 
measured in the units of Go = e'^/h and To = ^9 Go, respectively. The joint 
dimensionless cumulants 

Ke,m = {{iG/GoY{y/?or)) (2.3) 
§ The two can readily be extended to other symmetry classes, see Savin and Sommers (2006). 
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can be extracted from the expansion 

oo £ m 

logJ„{z,w)= J2 {~iy+"^^Ke,m, (2.4) 

v.rw. 

e,m=Q 

where kq.o = 0. In both Eqs. (2.2) and (2.4), the subscript n stands for 
n = min(iVL, TVr), and N = Nj, + is the total number of propagating modes 
(channels) in the leads. The notation § € CUE(iV) indicates that averaging runs over 
scattering matrices S drawn from the Dyson circular unitary ensemble (Bliimel and 
Smilansky 1990, Mello and Barangcr 1999, Mchta 2004). The latter is microscopically 
justified (Lewenkopf and Weidenmiiller 1991, Brouwer 1995) in the universal transport 
regime we are confined to. 

To perform the averaging in Eq. (2.2) in a most economic way, we employ a polar 
decomposition (Hua 1963, Baranger and Mello 1994, Forrester 2006) of S-matrix. 
Bringing into play a set of n transmission eigenvalues T = (Ti,---,T„) e (0,1)" 
distributed in accordance with the joint probability density function 

n 

Pr,{T)=c-'Al{T)l[T^, (2.5) 
this decomposition highlights Landauer's idea of viewing conductance as transmission, 

n 

G{T)=GoJ2t,. (2.6) 
Simultaneously, it reduces the expression for noise power [Eq. (1.4)] down to 

(n n \ 

J2^j+ffsY.Tj{l-Tj)\ . (2.7) 
i=i i=i / 

The parameter v in Eq. (2.5) is a measure of asymmetry between the leads, 
v = |7Vl — iVR,|, the notation A„(T) stands for the Vandcrmonde determinant 
A„(r) = YijKkC^k — Tj), whilst c„ is a normalisation constant. As the result, we 
are left with the JCGF in the form 

« n 

j„(z, w) = c- W n dT, t; r,,^(r,) (r), (2.8) 

■^(0,1)" jtTi 

where 

r,,^(T) = exp [-{z + w)T-wf0 T(l - T)] . (2.9) 

Although the above matrix integral representation of the JCGF 3^n{z,w) is by far 
more complicated than the one appearing in the integrable theory of conductance 
fiuctuations (Osipov and Kanzieper 2008), 

J„(^,0) = (exp(-^G/Go)>g^cuE(iv)=c;;' / H dT,- t; e"^^^ A^ (T), (2.10) 

it can still be treated nonperturbatively, much in line with the formalism used in the 

exact approach to zero-dimensional replica sigma models (Kanzieper 2002, Kanzieper 
2005, Osipov and Kanzieper 2007, Kanzieper 2009). 
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2.2. The T function theory of the joint cumulant generating function 

The "deform-and-study" approach (Morozov 1994, Adler et al 1995, Adler and van 
Moerbeke 2001) borrowed from the theory of integrable systems is central to the 
nonperturbativc calculation of 5'„(z,w). In the present context, the main idea of the 
method consists of "embedding" 5'„(z, w) into a more general theory of the r function 

1 r " 

Tr,{t;z,w) = - / \{dT^T^T,,^{T^)e''^*-'^^)^l{T) (2.11) 

-'(0,1)" 

which possesses the infinite-dimensional parameter space t = (ti, ^2, • ■ •) arising as the 
result of the t deformation 

oo 

V{t-T)=Y,tkT^- (2-12) 

k=l 

Studying an evolution of the r function in the extended {n,t,z,w) space allows us 
to identify various nonlinear differential hierarchical relations. A projection of these 
relations onto the hyperplane t = 0, 

n\ 

Jn{z,w) = —Tn{t;Z,w) , (2.13) 
c„ t=0 

generates, among others, a closed nonlinear differential equation for the JCGF 
J„(z,w). It is this equation [Eq. (2.24)] that will further supply the cumulants of 
noise power. 

The two key ingredients of the exact theory of r functions are (i) the bilinear 
identity (Date et al 1983) and (ii) the (linear) Virasoro constraints (Mironov and 
Morozov 1990). 

2.2.1. Bilinear identity and Kadomtsev-Petviashvili equation 

The bilinear identity encodes an infinite set of hierarchically structured nonlinear 
differential equations in the variables t = {ti,t2,...). For the model introduced in 
Eq. (2.11), the bilinear identity reads (Adler et al 1995, Tu et al 1996): 

/ rf.e«^(-*'-)r„(t-[z-])^II^±^±i^^ 

Coo 

= £ rf,e(°-^)-(^-*'^-)r^(f- - [z-^]) ^"+^J„*jjf (2.14) 

Here, a £ R is a free parameter; the integration contour Coo encompasses the point 
z = oo; the notation t ± [z~^] stands for the infinite set of parameters {tj ± z~^ /j}; 
for brevity, both z and w were dropped from the arguments of r functions. 

Being expanded in terms oit'—t and a, Eq. (2.14) generates a variety of integrable 
hierarchies (Osipov and Kanzieper 2009). One of them is the Kadomtsev-Petviashvili 
(KP) hierarchy. Its first nontrivial member, 

4 a2 fi2 \ /pa ^ 2 



is of primary importance since its projection onto t = [Eq. (2.13)] gives rise to a 
nonlinear differential equation for the JCGF 5'„(z, w). The resulting equation will 
further be used to determine the noise power cumulants we are aimed at. 
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2.2.2. Virasoro constraints 

Since we are interested in deriving a differential equation for 3^n{z,w) in terms of 
the derivatives over variables z and w, we have to seek an additional block of the 
theory that would make a link between tj -derivatives in Eq. (2.15) taken at t = and 
the derivatives over w and z. This missing block is the Virasoro constraints which 
reflect the invariance of the r function [Eq. (2.11)] under a change of the integration 
variables. 

In the present context, it is useful to demand the invariance under the set of 
transformations 

Tj -^fj+e ff\fj - 1), q> 0. (2.16) 

Employing by now a standard procedure (Mironov and Morozov 1990, Adler and van 
Mocrbckc 1995), one readily checks that the transformation (2.16) induces Virasoro 
constraints in the form 

[Lq+i{t)-L^{t)\Tn{t;z,w)=Q, q>0, (2.17) 

where a set of differential operators 

4(t) = t,{t) + 2f0 -[z + {l + M to]^ + (2.18) 
involves the Virasoro operators 

satisfying the Virasoro algebra 

[Lp, tq] = {p- q)tp+g, p,q>-l. (2.20) 
In Eqs. (2.18) and (2.19), the convention d/dto = n is assumed. 

2.2.3. Nonlinear differential equation for 3'n{z,w) 

To derive a differential equation for the JCGF 3^„(z,w), one has to project the first 
KP equation Eq. (2.15) onto the hyperplane t = 0. Spotting the identities {fjs > 0) 

d d , . 

—Tn[t;z,w) = - —Tn{t;z,w), (2.21) 

d d d 

f0-g^Tn{t;z,w) = —Tn{t;z,w)-{l + fp)—Tn{t\z,w), (2.22) 

we combine Eqs. (2.15) with the Virasoro constraints Eq. (2.17) taken at g = 0, 



^ . f d <9\„, { d d 
- + (1 + /.) H - ^) + (A^L + A^K) - iV.iVH 



Tn{t;z,w)=Q, (2.23) 
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to derive: 



q4 o2 

y^fffT^ + [2(^^L + Nn)U -2z + wil- /!)] ^ + 2(z ■ 



2m;) 



dzdw 



+3w 



d_ 

dw 



d_ 

dz 



d 



log9niz,w) + 6wffj[ — log3^n{z,w)] =0. (2.24) 



Owing to Eq. (2.4), this nonlinear equation considered together with the equation 
for 3^n{z,0) contains all the information about joint cumulants of the Landauer 
conductance and the noise power. The latter equation, written in terms of 



d 

cTn{z) = Nj^Nn + z—log3'n{z,0), 



(2.25) 



reads (Osipov and Kanzieper 2008): 



Q3 



& 



dz"^ 



d 



dz 



^o„—an{z) 



(2-(7Vl + A^r))'-47VlA^r 



— a„(^)-(A/-L + iVR-0)(T„(^)=O. (2.26) 

This can be recognised as the Chazy form (Chazy 1911, Cosgrove and Scoufis 1993) 
of the fifth Painleve transcendent 



52 



d 



dz 



d 



a„{z) + 2 ^a„{z) + (A/^l + ATr - ^)^a„(z) 



dz 



+4 



o-„(z) 



d 



A^L + ^a„(z) 



d_ 

dz 



On{z) =0 



(2.27) 



written in the Jimbo-MIwa-Okamoto form (Jimbo et al 1980, Okamoto 1987). For 
completeness, we have included a detailed derivation of Eq. (2.26) into Appendix A. 



2. 3. Recurrence solution for joint cumulants 

Combined with the cumulant expansion Eq. (2.4), the differential equation Eq. (2.24) 
furnishes the nonlinear recurrence for the joint dimensionless cumulants of 
conductance and noise power {(., m > 0): 

/| K.l+i,m-l + (1 - fp) l^e+2,m-l — 2 {Nj, + ATr) fff Ke+2,m 

-2 {£ + 2m + 1) Ke+i,m + (2^ + 3m + 2) Kt^m+i 

m — 1 



m 



i=0 



m — 1 
i 



E 



l^j+2,i Ke-j+2,m-i-l — 0. 



(2.28) 



To resolve it, one must know the boundary conditions whose role is played by 
cumulants Ki^o — {{{G/Go Y)) of the dimensionless Landauer conductance. These have 
been nonperturbatively calculated in our previous publication (Osipov and Kanzieper 
2008). Indeed, given the mean conductance 

Klfl = —r (2-29) 



and its variance 



(iVL + N^y 



(2.30) 
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the higher order cumulants «;^,o's are determined by the one-dimensional recurrence || 

[(TVl + Nnf - f] {£ + l)Ke+i,o + (A^L + Nn){2£ - 1) £Ke,o 

^ /£\ 
+£{£ -l){£- 2) K^_i,o - 2 ^(3j + l)(j -£?[■) «j+i,o «^-,-,o = 0. (2.31) 

Equations (2.28) and (2.31) represent the main result of our study ^. They 

provide a nonpcrturbativc description of the noise power fluctuations in the crossover 
region between the thermal and the shot noise (Savin, Sommers and Wieczorek 
2008) regimes (as discussed in the Introduction) by relating the temperature (9) and 
bias-voltage (v) dependent cumulants of the noise power to those of the Landauer 
conductance. 

Undoubtedly, the very existence of the above nontrivial relation (which 
emphasises a fundamental role played by Landauer conductance in transport 
problems) must be well rooted in the mathematical formalism and also have a good 
physics reason. As far as the former point is concerned, we wish to stress that a 
naive attempt to build a theory for the generating function ?'„(0,w) of solely noise 
power cumulants faces an unsurmountable obstacle: the KP equation [Eq. (2.15)] and 
appropriate Virasoro constraints [Eq. (2.17) at z = 0] cannot be resolved jointly in 
the hyperplane t = 0. This justifies the starting point [Eq. (2.2)] of our analysis. The 
physics arguments behind the peculiar structure of our solution are yet to be found. 



2.4- Noise power cumulants in the crossover regime 

Some computational effort is needed to read off explicit formulae for the noise power 
cumulants from Eqs. (2.28) and (2.31). Below we provide expressions for two families 
of joint cumulants expressed in terms of dimensionless cumulants k^,o = (((G/Go)^)) 
of the Landauer conductance. 

• Mean noise power 

P>)) = 46 Go [Nfp K2,o + Ki,o] (2.32) 

II Equations (2.29), (2.30) and (2.31) follow from the cumulant expansion 

log?-„(z,0) = ^(-l)^^«f,o, 
e=i ^■ 

see Eqs. (2.4) and (2.10), substituted into Eq. (2.27). Successive iterations of Eq. (2.31) yield the 
cumulants K£ q of Landauer conductance in the form 

(e - 1)! , , 

where the first few polynomials Pi{K) are: 

Pi{k) = K, 
P2{k) = K^, 



P3{k) = 4k^ - Nk^, 

P4{k) = 12(2- ^ - WNk^ + (AT^ + 1)«;2 



Here, N = N^ + N^. 

^ Notice that Eq. (15) in the paper by Osipov and Kanzieper (2008) contains typos. The correct 
formula is given by Eq. (2.31) 
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is generated by the lowest order member {£,m) = (0,0) of the recurrence 
Eq. (2.28). Being in concert with the known expression Eqs. (1.5) and (1.6), 
this result is a particular case of a more general formula 

J" 

Ke,i = Ke+ifi + ^Y+i (2.33) 



• Noise power variance 

= iA9 Go) 



2 \ f 2 / f 2 \ g 

-N^ -IJY '^4.0 + Nfi3 K3,0 + ( 1 + ^ j K2,0 " g /J/^a.o 

(2.34) 

is supplied by the (£, m) = (0,1) member of the recurrence. Its generalisation 
reads: 



2N^ \ f f I J \ 



j=0 

Here and above, A'' = A^l + A^r 



^^2l^ i' ( j) '*^J+2.o'^^+2-J.o- (2.35) 



The noise power cumulants ((y^)) of higher order {I > 3) can be calculated in the same 
manner albeit explicit expressions become increasingly cumbersome. Varying therein 
the parameters {9,v) from {0,0) to (0,w), one observes a smooth crossover between 
the thermal and the shot noise regime. 

2.5. Large-n analysis of joint cumulants: Symmetric leads 

The nonpeturbative solution Eq. (2.28) and (2.31) has a drawback: it does not provide 
much desired explicit dependence of conductance and/or noise power cumulants Ki^m 
on parameters of the scattering system. To probe such a dependence, we turn to the 
large-n limit of the recurrence Eq. (2.28). In what follows, the asymmetry parameter 
u will be set to zero. 

Under the latter assumption {v = 0), the joint cumulants Ki^m are solutions to 
the recurrence equation {£, m > 0) 



m 



fg K,e+4,m-l + (1 — ./I) Ki+2,m-l 



m+1 



-2 {£ + 2m+ 1) Ki+i^-m + {21 + 3m + 2) nt 

+ 6m /J ^ f . j X! ( • ) «^-j+2,m-i-i = (2.36) 

which must be supplemented by yet another recurrence {£ > 2) 
(4n^ - ^2) {£ + 1)k£+i,o + 2n(2£ - 1) £Ke,o 

+£{£ -1){£- 2) K^_i,o - 2 J2i^J + l)(i - ( • ) '^J+i.o '^(-jfi = (2.37) 

3=0 ^ 

that brings, in turn, a set of initial conditions K^fi to Eq. (2.36). 
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2.5.1. Cumulants of Landauer conductance 

It is instructive to start with the asymptotic analysis of Eq. (2.37). In case of sym- 
metric leads, the conductance cumulants of odd order vanish +, K2£+i,o = for all 
£ > 1 albeit Ki.o — n/2. As the result, one is left with the recurrence equation for the 

cumulants K2^,o of even order {£ > I) 

[4^2 - {2£ + 1)2] {£ + l)n2e+2,o + £{^£' - 1) «2£,0 

-8j2i^j + 2)(i - ef + 1 J '^2,+2,o K2e-2j,o = (2.38) 

subject to the initial condition [Eq. (2.30)] 

_ n"' _ 1 \ - 1 

- 4(4n2 - 1) ~ 16 ^^-^^^ 

Since, in the limit of a large number of propagating modes (n » 1), the conductance 
distribution is expected to roughly follow the Gaussian law (Politzcr 1989) with the 
mean k^ q = n/2 and the variance Kj.o = I/IQ [see Eqs. (2.29) and (2.30)], it is 
natural to seek a large-n solution to Eq. (2.38) in the form (j > 1) 

K2(.o = TF-^i'^ + '^'*2f.o, (2.40) 
Id 

where 5K2efi (with £ > 2) account for deviations from the Gaussian distribution. 
Putting forward the large-n ansatz 

o-=0 

where [see Eq. (2.39)] 

«2(2a) = (2.42) 
we further substitute it into Eq. (2.38) to derive: 

02^+2 (0) 1 (2£+l)' 



and 



Hence, 



02^0) 24(2^-1)!' 

04(0) _ 3! 
a2(2) ^ 



£>2, (2.43) 

(2.44) 



«2^(0) = -^y^- (2.45) 

Taken together with Eq. (2.42), this yields the leading term in the 1/n expansion 
[Eq. (2.41)] for conductance cumulants: 

l (2f-l )! 

4 (4n)2 

The higher order corrections to Eq. (2.46) can be obtained in a regular way. 



S^^^^^l-TJZylf- (2-46) 



+ At the formal level, this is direct consequence of the identity 3>i(z,0) = e "^?'„(— 2,0) holding as 
soon as = 0, see Eq. (2.10). 
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2.5.2. Dependence of the noise power cumulants on temperature and bias voltage 



Similarly to the previous subsection, we start an asymptotic analysis of the recur- 
rence Eq. (2.36) with singling out the large-n Gaussian part: 
n 



'm,0 



1 

16 



1 + ^ ) <5<;,o^m,i 



/I' 



(2.47) 



The Gaussian part was read off from Eqs. (2.33) and (2.35); the term Snt^m 
accommodates non-Gaussian corrections to the joint cumulants of Landauer 
conductance and the noise power. Their large-n behavior can be studied within the 

\/n ansatz 

1 at,m{o-) 



E 

o-=0 



(2.48) 



where a\fi{a) = 0. Substitution of Eqs. (2.47) and (2.48) into the two-dimensional 
recurrence Eq. (2.36) brings the recurrence equation {£-\-m> Q) 



m 



1 - 



(0) - 4/^a^+2,„(0) 



-2{e + 1 + 2m)a£+i,™(0) + {2i + 2 + 3m) a^,™+i(0) = (2.49) 

for the expansion coefficients a£,m(0) appearing in Eq. (2.48). The (unique) solution 
of Eq. (2.49), subject to the boundary condition 

«^,o(0) = ^^[l + (-l)^] 
derived in Sec. 2.5.1 [see Eq. (2.45)], reads {e + m> 0): 



(2.50) 



{£ + m-l)\ 



1 +(-1)^ 



//3 



- 1 



(2.51) 



2 ' J ' ' ' \ 2 
Combined with Eq. (2.48), it yields the leading term in the 1/n expansion for joint 
cumulants of Landauer conductance and the noise power: 



m 



1)! 



k 

2 



1 



2 



1 



(2.52) 



^ (4n)^+™ 

Equations (2.47) and (2.52) are the central result of this Section. 

In particular, it brings an the following large-n expression for the cumulants of 

noise power in case of symmetric leads: 



i^)) (Goer 



2n 



1 + ^ + I 1 + 



+ 



(^-1)! 
8n^ 



ll - I 
2 



(2.53) 



Explicit dependence of the noise power cumulants on both the temperature 9 and 
the bias voltage v enters through a single function [see Eq. (1-7)] that depends 
on the ratio /3 = v/26. Based on Eq. (2.53), it can further be shown that small 
but nonvanishing cumulants of the third and higher order arc responsible for long 
exponential tails in the otherwise Gaussian distribution of the noise power (compare 
with Vivo, Majumdar and Bohigas 2008). 
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3. Conclusions 

In summary, we have presented an advanced formulation of the recently proposed 
integrable theory of quantum transport (Osipov and Kanzieper 2008) to study 
statistics of noise power fluctuations in a chaotic cavity with broken timc-revcrsal 
symmetry in the crossover regime {9, 0) (0, v) between thermal and shot noise. 
By relating the cumulants of noise power to those of the Landauer conductance, we 
determined dependence of the noise power cumulants (as well as of joint cumulants 
of Landauer conductance and the noise power) on the bias voltage temperature 0, 
and the number of channels A^l.r in the leads attached to a cavity through ballistic 
point contacts. 

We are confident that ideas of integrability combined with the scattering matrix 
approach are able to provide a nonperturbative description of many more transport 
phenomena in chaotic cavities. Quantum transport in cavities with losses (Doron, 
Smilansky and Frenkel 1991, Beenakker and Brouwer 2001, Simon and Moustakas 
2006) and non-ideal leads (Brouwer 1995) are just two examples of chaotic scattering 
systems whose detailed study is much called for. 

Note added. — Recently, we learnt about the paper by Khoruzhenko et al (2009) where 
an alternative approach was developed to describe statistics of conductance and shot- 
noise power in chaotic cavities with and without time-reversal symmetry. Possibly 
triggered by the earlier paper by Novaes (2008), these authors combine a theory of 
the Selberg integral with the theory of symmetric functions to evaluate the (joint) 
moments of Landauer conductance and the shot-noise power in terms of series over 
all partitions of the moment's order. In particular, Khoruzhenko et al (2009) confirm 
our large-n formulae Eqs. (2.46) and (2.53) at zero temperature. 
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Appendix 

A. Cumulant generating function for the Landauer conductance and the 
fifth Painleve transcendent 

Relation to a gap formation probability. — The "simplest" (albeit not operative) way to 
observe that the conductance cumulant generating function 5'„(z, 0) [Eq. (2.10)] can be 
expressed in terms of Painleve V is to spot that 5F„ (z, 0) is essentially the gap formation 
probability within the interval (z, +oo) in the spectrum of an n x n Lagucrre unitary 
ensemble. Indeed, the transformation of integration variables Xj = zTj in Eq. (2.10) 
yields 




n 



(A.1) 
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A nonperturbative evaluation of the above n-fold integral is readily available (Tracy 
and Widoin 1994, Forrester and Witte 2002) eventually resulting in the following 
Painleve V representation (Osipov and Kanzieper 2008): 

,o'n(i) - n{n + i/)' 



J„(^,0) = exp 



dt- 



t 



(A.2) 



Here, an{t) satisfies the Jimbo-Miwa-Okamoto form of the Painleve V equation (Jimbo 

et al 1980, Okamoto 1987): 

(to' + [cJn - ta'„ + 2{a'„f + (2n + ,^)a'„]^ + 4{a'J^a'„ + n){a'„ + n + ly) = (A.3) 

subject to the boundary condition cryit — > 0) ~ n{n + z/). Keeping in mind the pa- 
rameterisation n = m.m{Ni,, Ni,) and v = \Ni, — JVrI, one concludes that Eq. (A.3) is 
equivalent to Eq. (2.27) announced in Sec. 2.2.3. 



Direct evaluation of 3^n{z,Q)- — To directly evaluate 5'„(2;,0) defined by Eq. (2.10), 
we introduce the associated r function 



1 r " 

Tn{t; z) = -l n e-^^+^(*=^^) A^(T) 



such that 



J„(^,0) = —Tn{t;z) 



t=0 



and make use of the KP equation Eq. (2.15), 



dtf^^dq 



dtidts 

supplemented by the Virasoro constraints * 

[Lg+i{t)-Lg{t)]Tnit;z) = 0, q>0 
where a set of differential operators 



0, 



d 



dta 



_d_ 



(A.4) 

(A.5) 

(A.6) 
(A.7) 
(A.8) 



iq+l ^uq 

involves the Virasoro operators Eq. (2.19) [the convention d/dto = n is assumed]. 

In order to project the KP equation Eq. (A.6) onto t = 0, we need only two 
Virasoro constraints labeled by g = 0, 



and q = 1 

oo 



d 



dti 



d_ 



9 ,^ , d 

.- + (2n + . + .)- 



logr„(t; z) = n(n + i/), (A. 9) 



dti 



■j+2 



dtj+i 



-{2n + u) 



9 ,^ , d 

- z— + {2n + u + z) 

Ots 



dt2 



_d_ &^ 
dti ^ dti 



logr„(t;^) 



dti 



log T„(t,^) =0. (A.IO) 



* Equations (A.7) and (A.8) readily follow from Eqs. (2.17) and (2.18) upon setting = and 
w = 0. 



References 



14 



Here 

d d 

— r„(t; z) = -TrT„(t; z). (A.ll) 
oti oz 

Lengthy but straightforward manipulations with Eqs. (A. 9) and (A. 10) as well as with 
their derivatives over ti and t2 projected onto t = result in Eqs. (2.25) and (2.26). 
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